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Abstract 

We have optimized the zero frequency first hyperpolarizability (3 of a one-dimensional piecewise linear 
potential well containing a single electron by adjusting the shape of that potential. With increasing num- 
bers of parameters in the potential, the maximized hyperpolarizability converges quickly to 0.708951 of 
the proven upper bound. The Hessian of ft at the maximum has in each case only two large eigenvalues; 
the other eigenvalues diminish seemingly exponentially quickly demonstrating a very wide range of nearby 
nearly optimal potentials, and that there are only two important parameters for optimizing j3. The shape 
of the optimized wavefunctions converges with more parameters while the associated potentials remain 
substantially different, suggesting that the ground state wavefunction provides a superior physical descrip- 
tion to the potential for the conditions that optimize the hyperpolarizability. Prospects for characterizing 
the two important parameters for near-optimum potentials are discussed. 



I. INTRODUCTION 



The non-linear response of a electron to an electric field is an important fundamental and 
practical issue in non-linear optics [HE]- In a notable series of papers, Kuzyk et al. have shownp]- 
5j that all of the measured off-resonant first hyperpolarizabilities (5 ZZZ — specifically the second 
derivative of the polarization in a given direction z with respect to an electric field with vanishingly 
small frequency and wavevector in that direction — of known molecules are considerably smaller 
than a theoretical upper limit f3 max derived from the Thomas-Kuhn sum rules and the sum-over- 
states formula 



Here N is the number of electrons participating and E w = E\ — Eq is the energy difference between 
the ground state and first excited state. This upper bound need not be realizable as it is unclear 
that the matrix elements required can be achieved. 

One possible strategy to achieve large hyperpolarizabilities is to optimize the shape of the 
potential in which the electrons are confined. Kuzyk et al. have studied a sub-problem of this, 
namely performing the optimization for the single-electron case with both one and two dimensional 
potentials. They suggest that implementing the qualitative (modulated) features observed in their 
optimized potentials might allow design of molecules with higher hyperpolarizabilities|U [5]. While 
these calculations have resulted in a number of potentials with large hyperpolarizabilities, they do 
not distinguish which features of the optimized potential are required for high (5 and which are 
artifacts of the minimization procedure or chosen parametrization. 

In this work we have re-examined the optimization of the first hyperpolarizability for a single 
electron moving in a potential in one dimension. We find convincing evidence that the best possible 
hyperpolarizability is substantially smaller than previously published limits: instead we propose 
a new lower actual maximum of approximately 0.708951/3 max . Moreover, we find a surprisingly 
broad class of potentials that have values of essentially indistinguishable from the maximum 
and that, in effect, there are only two important constraints on potentials to have near-optimal 
P. Perturbations to an optimized potential that do not affect its hyperpolarizability include those 
where the ground state wavefunction is small as well as rapid changes where the ground state 
wavef unction is large. 

To quantify the effect of such perturbations, we present the problem in the language of opti- 
mization theory. Here, one speaks of maximizing an objective function — in this case the "intrinsic" 
hyperpolarizability (3 int = f3/(3 max — with respect to some set of parameters, i.e. those that spec- 
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ify a potential. Once an optimum value has been found, it is natural to ask how sensitive that 
value is to variations in the parameters. This information is contained in the Hessian matrix, the 
eigenvalues of which are essentially the local principal curvatures of the objective function at the 
maximum [7J. 

It is quite possible, and will be shown to be the case in the present work, for these eigenvalues to 
have very different magnitudes. While the eigenvalues of the Hessian lack a physical interpretation, 
their importance for the numerical optimization procedure is immediate: large ratios of eigenvalues 
of the Hessian mean that the maximization procedure is, in effect, attempting to find the highest 
point on a very narrow ridge with a flat top. Generally, optimization programs will stop when 
improvement of the function is slow, or the approximate gradients are small. Alternatively, the 
program may stop after a specified number of iterations; the result may be accepted if the value of 
the objective function has converged even if the values of the parameters have not. It is therefore to 
be expected that the numerical errors in parameter combinations associated with small eigenvalues 
of the Hessian are much larger (by some inverse power of the eigenvalue roughly between —1/2 and 
— 1) than those associated with large eigenvalues [7J. While these estimated errors, unfortunately, 
are typically not reported by such programs, examining them allows one to be confident both 
that the actual maximum hyperpolarizability has been found and that the potentials are not 
inappropriately over-specified. For the present work, they also provide information on the likely 
errors in reported potentials. Since very small eigenvalues indicate a flat maximum, small errors 
associated with discretization, numerics, etc may substantially change the result. 

There are multiple possible causes for large ratios of eigenvalues. It may be that some funda- 
mental feature of the problem is responsible, or that the representation chosen is deficient in that 
it admits this behavior. In the present case, we expect that the hyperpolarizability depends on 
the potential primarily where the ground state wave function is substantial, so parameters that 
control the shape of the potential elsewhere are likely to affect the hyperpolarizability negligibly. 
Moreover, since the hyperpolarizability has been shown to depend only on dipole matrix elements 
and may be derived entirely as a multiple integral of the ground state wavefunction [8j, changes 
in the potential where the ground state wavefunction is large are also irrelevant if they tend not 
to alter the wavefunction. For example, rapid variations in the potential which approximately av- 
erage to zero will not significantly affect the hyperpolarizability. For numerical optimizations like 
those performed below, it is desirable to specify the potential in such a way as to avoid including 
parameters made redundant by these two effects: the fewer irrelevant parameters included in such 
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optimizations, the greater the precision that can be achieved for the remainder. 

The number of parameters that characterize the chosen potential is also important: if too few 
coefficients are included the problem is over-constrained and the optimal value of /3j n t will be 
significantly less than the true maximum. Conversely, if too many coefficients are included the 
problem might become imder-constrained. In such a case many potentials yield comparably opti- 
mal values and interpretation of the potential requires care: important features may be obscured 
by irrelevant ones associated with small eigenvalues of the Hessian. Moreover, the convergence 
will be slow particularly for numeric schemes, like those used to date on this problem, that only 
evaluate the objective function and not its derivatives. 

A known technique[9j for investigating such optimization problems is to increase the number of 
parameters describing the problem and examine the eigenvalues and eigenvectors of the Hessian 
at the maximum with each additional parameter. To investigate whether the potential is under- 
constrained we have chosen the parametrization as suggested above and varied the number of 
parameters that characterize the potential. In each case, we optimized /3j nt and calculated the 
eigenvalues and eigenvectors of the Hessian at the maximum. This analysis allows us to state 
which aspects of the potential are largely irrelevant to the hyperpolarizability. It also reveals 
that effectively only two continuous parameters need to be specified to determine near-optimal 
potentials in the present problem. 



II. MODEL 



There are myriad possible parametrizations for ID potentials and several quite simple potentials 
have hyperpolarizabilities that are a large fraction of the putative maximum, e.g. the semi-infinite 
triangular well and the clipped harmonic oscillator [5J. We have chosen to represent the potential 
as a piecewise linear function that consists of iV + 1 segments 



V(x) 



AqX + B X < x 

A n x + B n x n ^i < x < x n , n e {1, 2, N - 1} 
A N x + B N x>x N _i 



(2) 



with positions x n and slopes A n as the parameters. The x n are strictly ascending and this must 
be enforced in the optimization. Since /3 is invariant under both translation of x and the addition 
of a constant to V(x), two constants are fixed Xq = and B = with no loss of generality. This 
specificity is desirable: it decreases the dimensionality of the problem that must be solved in order 



to find a given optimum potential, eliminates part of the null space of the Hessian and so increases 
the typical size of the gradient during an optimization. Furthermore, it is substantially easier to 
compare optima and verify that the global optimum is unique. In order to enforce continuity, we 
take B\ = B$ and the remaining constants B n are given by 

n-l 

B " = ^2( A m+i ~ A m )x m , n > 1. (3) 

m=l 

A remaining degree of freedom irrelevant to the determination of (3 max , the energy scale associated 
with the potential, is removed by setting A\ = 1. This does restrict the class of potentials repre- 
sented to those in which the leftmost linear segment has negative slope and the second has positive 
slope, and so for iV > 1 the effect of choosing A\ = — 1 was investigated and no improvement in 
the optimal hyperpolarizability was found. In order to ensure a bounded wavefunction it is also 
necessary to impose A < and A N > 0. Having imposed these constraints, there remain 2N — 1 
free parameters Pj = {Aq, A 2 , A^, Xi, x^-i}- In order to allow the possibility of having an 
even number of free parameters, the optimizations were performed in increasing order of parameter 
number and for even cases the left hand slope A was fixed to the value obtained from optimizing 
the potential with one fewer parameter. 

The hyperpolarizability for a given potential was calculated by the following procedure: The 
Schrodinger equation for one electron in the n-th segment is, in units where e = 1, h = 1 and 
m e = 1, 

r 1 A 2 1 

i>n = Eij n (4) 



Id 2 

2dx 



2 + (A n + e)x + B n 



where an electric field of strength e is also included since it is desired to calculate {3 at e = 0. Its 
solutions are the well-known Airy functions 



ip n (x) = C n Ai 



^2(B n -E + x(A n + e)) 



+ £„Bi 



V2(B n - E + x(A n + e)) 



(5) 



(A n + e) 2 /3 

The 2N + 2 constants C n and D n are determined by imposing continuity of tp(x) and its derivative 
ip'(x) at each of the set of points x i: yielding iV pairs of equations. The two remaining equations 
are found by requiring the wavefunction to vanish as x — > ±oo, mandating that D = = 0. 
For the case where the number of free parameters N = 2, the potential is infinite for x < Xq and 
the two equations at x = x are replaced with the usual requirement that the wavefunction vanish 
there. All of these equations are linear in the parameters C n and D n and may be represented in 
matrix form 

W ■ u = (6) 



where u is the vector of coefficients {C n , D„} and the matrix W depends on E and e. Solutions to 
(g exist if 

det W = 0, (7) 

which was readily solved numerically with e = to yield an ordered set of energy levels {E{\. The 
hyperpolarizability was then calculated as follows: The Jacobi formula for the derivative of the 
determinant is 

d detW/ = Tr (adjW- ^] (8) 



de \ de 

where adjiy is the adjugate of W . Since det W = 0, for all e, 4 det W = 0. Using the chain rule 

dW _ dW dW dE 

~de~ ~ ~de~ + d£ de~' ( ' 

the derivative of the energy with respect to e may be readily calculated 

dE _ Tr(adjiy-f^) 

de Tr(adj^-ft)' 1 Uj 

Higher derivatives may be obtained by differentiating ([8]); performing this once yields an expression 



d 2 E Tr [| (adjW) • ^ + adjlU ■ W] 



de 2 Tr(adjW/-f|) 

where 

<9 2 iy <9 2 iy 



de 2 de 



2 dEde + de dE 2 



(12) 



and a second time yields ^-f 



d 3 E _ Tr(^ (adjWO " + 2| (adjW) ■ § + adjW/ • 
de 3 " " Tr(adjH/- H) 

where 



(13) 



_ 9 3 H/ d 2 £/<9 2 W/ d£9 2 lf\ d£ / d 3 W 
~lk 3 ~ + de 2 " \dEde + ~de~~dE 2 ) + de~ \dEde 2 + 



dE d 3 W \ fdE\ 3 d 3 W . iN 

+ t ( 14 ) 



de d£ 2 de ) \de ) dE 3 



and then the hyperpolarizability /3 is 



(15) 

E=E ,e=0 



2 de 3 

from which the intrinsic hyperpolarizability 

Pint = P / Pmax (16) 

can be calculated using Q together with the numerically determined E and E\ from 
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III. RESULTS AND DISCUSSION 



A program was written in Mathematica 7 using the above formulation to calculate the intrinsic 
hyperpolarizability for a potential with a variable number of parameters and with arbitrary values 



for those parameters, expressions for the necessary derivatives in (13) being evaluated automat- 
ically. The program is provided as supplementary material to this work. For a given number of 
parameters ranging from two to seven, the generic maximization routine FindMaximum with the 
Interior Point method and numerically calculated gradients was used to find locally optimal values 
of Pi n f, in each case the maximization was begun from many different randomly selected starting 
points in the parameter space and the best obtained locally optimized values of (3i n t taken to be 
the global maximum. We found for larger numbers of parameters that it was necessary to increase 
WorkingPrecision to beyond machine precision. A unique global maximum, to numeric preci- 
sion, was obtained for each parametrization independent of the starting point of the optimization. 
Care was taken that the optimization had converged to a numerically accurate result for each 
parametrization. At each of these globally optimal values, the Hessian matrix of the objective 
function 

H ti =^( m (it) 



dPidPj VA, 

was calculated numerically with second-order finite differences and its eigenvalues hj and associated 
eigenvectors v\ found. For each given number of parameters, the optimized 0i n t are displayed in 
table [I] together with the ratio of the smallest to largest eigenvalue of the associated Hessian. 
The optimized potentials are plotted in fig. [T] together with the wavefunctions for the ground 
state and first excited state. Here the wavefunctions have been normalized and the potentials and 
wavefunctions have been shifted and rescaled so that (x) = and E\ — E = 1, using 



x — (x 



- (xMEk - E Q f/\ V(x, {P}) = (V(x, {P}) - E )/(E 1 - E 



This allows us to compare transparently, without the confusion of irrelevant scaling, the wave- 
functions and potentials for different parametrizations of the potential in order to try to understand 
which features are pertinent to optimizing f3 int . Note that the optimal potentials for different num- 
bers of parameters are markedly different despite the fact that additional parameters generalize 
the potential. In contrast, changes in the wavefunctions as the number of parameters increases 
are far more subtle. As anticipated, the adjustable features in the optimized potentials are found 
to be strictly in the region where the ground and first excited state wavefunctions are large, even 
though this was not enforced in the optimization. 



Like previous analyses, we find the very best value of f3i nt to be somewhat lower than 1 which 
corresponds to a /3 equal to the theoretical upper limit of ([I]). With only four parameters in the 
potential, we find a value of (3i nt = 0.708764 comparable to the best previously known, and the 
highest value so far obtained (3i nt = 0.708951 with only 6 parameters — though it is a refinement 
of only 0.02% over previous results. The addition of a seventh free parameter yields no further 
increase in (3 int to within our numerical accuracy. We believe that this value is very close to the 
actual maximum achievable for a single electron moving in one dimension, or more speculatively, 
for the hyperpolarizability of a single electron with all components of the electric field along a 
single direction as e.g. f3 zzz . Formally it is simply a lower bound for the actual upper bound on 
Pint- Given that we have carefully verified that the optimization has converged, it is clear that 
there are no potentials similar to these that have f3 significantly higher by more than a fraction of 
a percent. It remains possible that there exists a potential with still larger hyperpolarizability that 
is not accessible with our parametrization, or which is qualitatively different from the potentials 
we have explored. However, to the best of our knowledge, all potentials that have been shown to 
have Pi n t approaching this value have ground state wavefunctions similar to those we have found. 
The wide range of parametrizations that have been explored, particularly by Kuzyk et al. [2H5], 
strongly suggests that qualitatively different but substantially better potentials do not exist and 
that the true optimum hyperpolarizability is close to those presented here. 

For our optimized potentials, there is a substantial variation in the eigenvalues of the associated 
Hessian as may be seen from table [I] and fig. [2] Surprisingly, it appears that effectively only 2 
parameters, those associated with the two large eigenvalues in the Hessian, need to be adjusted 
in order to achieve the maximum hyperpolarizability. There are, therefore, a very wide range of 
potentials near to the optimized one that have hyperpolarizabilities close to the maximum. This 
is consistent with the speculation of Kuzyk [2J that only some aspects of the first few wavefunc- 
tions — specifically the ratios of energy differences and dipole matrix elements of the ground and 
first two excited states — significantly affect the objective function: relatively few parameters in the 
potential are required to adjust these three low-lying wavefunctions accurately enough to specify 
quantities derived from their integral such as dipole matrix elements. Variations of the potential 
where these wavefunctions are small, or rapid variations where they are large, have little effect on 
either energies or dipole matrix elements. 

The eigenvectors of the Hessian allow us to determine what aspects of the optimized potential 
are pertinent to maximizing the hyperpolarizability. The variation in the potential AV^(x) in the 
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direction of the eigenvector associated with the j-th largest eigenvalue of the Hessian h? is 

dV(x, {P t + avj}) 



AVHx) 



da 



(19) 



where v\ is the i th component of the j-th eigenvector jTT] . A small variation in the potential 
along this direction V — > V+aAV^ will result in a corresponding change in the hyperpolarizability 
Pint ~^ Pint + \h? a 2 . The significance of such a variation is obscured by the fact that it typically 
contains some component that serves only to translate and rescale the energy and length. Since 
the hyperpolarizability is not affected at all by these operations, it is desirable to remove their 
associated components from the AVK This may be achieved by rescaling V and x in the right 



hand side of (19) as a function of a to preserve the properties (x) = and E\ — Eq = 1 using 



(18). The resulting adjusted AV J (x) no longer contain irrelevant rescalings or translations but 
appear to decrease in apparent size with decreasing eigenvalue (fig. [3]), indicating that they are 
increasingly associated with either the irrelevant quantities (x), E and E\ — E or with changes 
where the wavefunctions are small. Unfortunately, those AV^{x) associated with large eigenvalues 
do not seem to converge with increasing numbers of parameters and so reveal little about which 
features of the potential are truly important to the hyperpolarizability. 

This suggests our parametrization of the potential could be improved, and also that a physically 
more relevant measure of "distance" between potentials is desirable. This measure is implicit in 
the process of finding the eigenvalues and eigenvectors of the Hessian; in the above analysis, it had 
the form 



IIV^ — Vbll 2 = \Pi-P 2 \ 2 (20) 

where the Pj is the vector of parameters that specifies the i-th potential. While this idea of 
distance — referred to hereafter as the "numerically natural measure" — matters for the numerical 
process of finding the maximum, it is largely irrelevant to the design and synthesis of molecules 
that have large hyperpolarizabilities. To explore the possibility that the idea of distance between 
potentials matters for the interpretation of what matters for the potentials, we have constructed 
more physical measures of the distance between two potentials V\ and V2 as follows. First rescale 



each potential using (18) to set (x) = E = and E\ — E = 1, yielding scaled potentials Vf, V{ 
for which the quantities irrelevant to the hyperpolarizability have been adjusted to have standard 
values. The measure is then 
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oo 




dx Pk {x)\V°-V> 



(21) 



— oo 



where p k is a positive-definite weighting function. There are a wide variety of plausible possible 
choices for p k . we used po = ipQ and p^o = V'o + ip\ where tpo and ipi are the normalized, scaled 
wavefunctions for the optimized potential as these measures serve to remove physically irrelevant 
differences between potentials, and are consistent with a priori expectations that only changes in 
the potential where low lying wavefunctions are large matter to the hyperpolarizability. 

Eigenvalues of the Hessian were recalculated in each of these measures and are displayed in 
table [IIJ Clearly, the eigenvalues and eigenvectors of the Hessian are quite different for the different 
ways of measuring distances; unfortunately, this is particularly true for the eigenvectors that are 
associated with large eigenvalues of the hessian, i.e. those eigenvectors that tell us which changes 
to the potential do matter to the hyperpolarizability. Physically this is clear: changes that do 
not change the hyperpolarizability can be freely included in changes to eigenvectors with large 
eigenvalues if they do not also increase the distance. To examine this more precisely, consider the 
generalized eigenvalue problem with different measures of distances that may be written as 



Here, in a given parametrization, H k j is the hessian of the objective function (Pint)', M k j is the 
positive definite hessian of the measure of distance in the same parametrization; the v\ are the 
eigenvectors and the h? are the eigenvalues. As only the eigenvectors and ratios of the eigenvalues 
are of interest, M can be chosen to have unit determinant. Some idea of the importance of the 
measure can be obtained by considering two nearly identical measures and calculating the changes 
in the eigenvalues and eigenvectors using perturbation theory in 5M = Mi — M 2 . Let h? a and 
be the eigenvalues and eigenvectors within measure M a where eigenvectors are normalized by 
v^M a v^ = 5jk and 5M = (Mi) -1 / 2 5M (Mi) -1 ' 2 . Supposing the eigenvalues h{ and eigenvectors 




H ki vl = h j M ki vl 



(22) 





v{ for the measure Mi are known and that 51SAj k = vjiJMvf is small for j ^ k, the eigenvalues 
and eigenvectors h 3 2 and v 2 in the second measure M 2 may be found from a routine perturbation 



calculation: 




(23) 



and 




(24) 



10 



If h{ 3> h\ any matrix element £Mjfc between these eigenvectors will result in a comparably sized 
change in the eigenvector: in other words, the contribution of an eigenvector that corresponds to 
a small eigenvalue to an eigenvector with a large eigenvalue is controlled largely by the measure 
and not by the nature of the hyperpolarizability. If, on the other hand, h\ <C h{, changes in the 
measure result in rather little change to the eigenvector. Unfortunately, the eigenvectors for which 
h\ ^> h\ are the ones that ought to indicate which features of the optimized potential most affect 
the hyperpolarizability. For this reason, the question of which features of the potential do matter 
to Pi n t is much harder to answer than the matter of what does not. 

The changes in the scaled potential associated with each of the eigenvectors of the Hessian are 



displayed in fig. |3jwith the numerically natural measure of eq. 20 and the measure of eq. (21) 
with po = "00 5 a l so plotted are the corresponding changes in the ground and first excited state 
wavef unctions. Similar calculations for the other measures proposed produced similar behavior. 
Again, the superiority of a wavefunction description of the problem is clear from the fact that 
changes in the wavefunctions for the smaller eigenvalues are substantially smaller than those in 
the potential. 

The scaled h? and AW(x) within a physically- motivated measure unambiguously identify the 
range of potentials close to the maximum. The breadth of this range with a variety of measures of 
distance suggests that the task of synthesizing near optimum chromophores or, at least, designing 
near-optimum potentials, may be easier than has been previously supposed; on the other hand 
it appears that while the parametrizations of the potential used previously have yielded useful 
estimates of the upper limit on the hyperpolarizability, they have also tended to obscure the 
underlying physics. 

To further illustrate this, and to attempt to find a clear exposition of the important constraints 
for optimum hyperpolarizabilities, we calculated and display in table [I] two other physical quanti- 
ties: the dipole transition matrix element from the ground to first excited state x 10 = (1 |x| 0) and 
the change in dipole moment between the ground and first excited state, Axio — (1 \%\ 1) — (0 |x| 0). 
Both of these have been previously identified|l] from the sum-over-states formula as being impor- 
tant to optimizing j3. Importantly, they can also be measured. It is apparent that both quantities 
converge with increasing number of parameters, X\q to ~ 55% of the maximum possible value while 
less rapidly Axiq ~ 1.26. Thus there appears to be, consistent with earlier work, a rather specific 
optimum value of x 10 and a less clear optimum choice for Axiq. To further understand this we 
examined the derivatives of these parameters in the direction of each of the various eigenvectors 
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of the hessian. These are given for N = 7 parameters in table [IT] and clearly diminish with the 
associated eigenvalue. The largest eigenvector with the largest eigenvalue seems quite strongly 
to constrain a linear combination of them, which is largely Axio and the next largest eigenvalue 
constrains another combination that is largely Xiq. The fact that Axio is largely controlled by 
the largest eigenvalue is unusual given that X\o converges more quickly thanAxio- Nevertheless 
these quantities appear at first sight to offer a superior indication of how far a particular po- 
tential is from the optimum because deviations from the optimum alter them linearly, while the 
hyperpolarizability is only altered quadratically. It may be that combinations of these quantities, 
together with some other quantity presently unknown form a better indicator of the closeness of 
a potential to optimum, particularly for the second largest eigenvalue for the derivatives are not 
terribly large. Moreover, such ideas must be interpretted with care. It is, for example, physically 
clear that constructing a potential with a large barrier between the states that are important to 
the hyperpolarizability and a physically distant state precisely resonant with the first excited state 
will result in first and second excited states that are superpositions of the important excited state 
and this distant state. This will (of course) dramatically change x\o and Axio without signifi- 
cant effect on the hyperpolarizability. Nevertheless, for the single electron problem, and absent 
motions in special directions, these parameters, with the specific targets calculated herein seem 
reasonable heuristic, experimentally accessible measures of the distance of a physical system from 
the optimum. 



IV. CONCLUSIONS 



We have determined that only 2 parameters are important in optimizing the zero frequency 
non-resonant hyperpolarizability of a single electron in a one- dimensional potential; these are 
associated with large eigenvalues of the Hessian of the hyperpolarizability, which are dramatically 
larger than the other eigenvalues with any of a number of physically motivated measures. Even 
with our careful optimization with few, appropriately chosen, parameters it is very difficult to 
distinguish the truly optimal potential from potentials with hyperpolarizabilities very close to 
the maximum. This is clearly a generic problem: there many, many potentials with virtually 
identical hyperpolarizabilities that will appear to most observers quite different from the true 
optimum potential. It seems relatively easy, particularly in light of the wide variety of such 
potentials obtained by others[3H5], to design and produce systems whose potentials result in fli nt 
near the apparent maximum. However, the specific features of potentials that result in such 
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hyperpolarizabilities are, almost in consequence, difficult to derive and explain. Our calculation is 
sufficient to reveal a wide variety of things that do not matter to the optimization, and to make 
clear that the optimum is more clearly expressed in terms of the wavefunction than the potential. 
While this is a significant improvement over prior works that take no special steps to understand 
which features of the potential are irrelevant, we were still unable to characterize to our satisfaction 
in a concise manner what does matter about optimum potentials. We believe that a quite different 
calculation — described in part below — would be required to determine this, though not uniquely 
as the question of relevance depends in part on specifying a measure of distance between potentials 
and these are motivated by a particular chemical or physical realization. 

Given present results we are exploring other parametrizations of the potential: a judicious 
choice would only include variable features relevant to the maximization problem. The present 
work shows both that the optimized potentials become smoother with additional parameters and 
that rapid variations in the potential are irrelevant; this implies that the actual optimum potential 
is analytic. A promising approach is to expand the potential in a series of analytic functions such 
as Hermite polynomials and to maximize the hyperpolarizability with respect to the coefficients. 
Such a parametrization has some attractive features: If the optimal potential is analytic, the 
expansion should converge quickly with the higher order polynomials clearly corresponding to 
rapid variations in the potential that are known to be irrelevant. The following observations: 
that additional parameters make the potentials smoother; the eigenvalues of the hessian decrease 
rapidly; and that the eigenvectors somewhat resemble hermite polynomials with two added to the 
order; all suggest that the best potential is analytic as these characteristics would be expected 
for analytic potentials. Second, the numerically natural measure of distance between potentials 
induced by the weighting function of the Hermite polynomials is a physically reasonable one, 
since the weighting function itself resembles the electron's probability density in the ground state. 
We therefore believe this parametrization ought to reveal more rigorous lower bounds on the 
maximum possible hyperpolarizability and a clearer description of the optimal potential. It is also 
expected, for numerically feasible optimizations, to appreciably decrease errors in the eigenvalues 
and eigenvectors of the hessian from purely numerical sources such as discretization and truncation. 
This, in turn, may make physical interpretation of these quantities, and specifically the parameters 
that most matter to the hyperpolarizability appreciably easier. 

The consequences of our main results for the design of new non-linear chromophores are less 
obvious. Our clearest result is that numerically optimized potentials, such as those obtained in this 
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work and elsewhere[U [5], need very careful interpretation prior to use as design paradigms: both 
the parametrization and the optimization method are likely to result in appreciable "noise" that 
is without physical consequence. It seems remarkably difficult from our (or other) work to specify 
the range of nearly optimal potentials per se in any way that will matter to practical design. The 
N = 2 and N = 7 potentials of figure 1 differ by only slightly more than 1% in the optimum 
hyperpolarizability but appear physically dramatically different. 

We have also demonstrated convincingly that the true upper bound on the hyperpolarizability 
of a single electron is appreciably below the bound derived from the dipole matrix element sum 
rules. With n electrons, the maximum possible hyperpolarizability is predicted to be n 3 / 2 times 
that for a single electron. In fact this can be achieved in a system in which there are strong 
attractive interactions between the electrons that consequently move as a single quasi-particle in 
a potential similar that calculated above [10]. This mechanism is likely only to be realizable in 
exotic systems due to the Coulomb repulsion of the electrons. A more practical lower limit on 
the bound is n times the single electron maximum hyperpolarizability, which would correspond 
to the situation where the electrons are confined to non-interacting wells. Including repulsion, 
we speculate that the sum-rule bound increasingly overestimates the true bound with increasing 
numbers of electrons: two non-interacting electrons with opposite spins and Fermi statistics have 
twice, not 2 3 / 2 times the upper bound on the susceptibility derived from the sum rules, while n 
non-interacting bosons would have only n rather than n 3 ^ 2 times the one particle bound. The fact 
that the bound for non-interacting electrons is smaller than that for attractive electrons suggests 
that it may also be a bound on electrons that repel each-other. Hence, it would appear that extant 
materials fall less spectacularly short of actual bounds than is presently claimed, although there 
is likely still substantial improvement possible in the hyperpolarizabilities. 

Most importantly, it appears that surprisingly few parameters of the potential need to be ad- 
justed to achieve close to optimum hyperpolarizability. Thus, while we do believe that significantly 
better chromophores can still in principle be made, it seems likely that careful, multidimensional 
optimization of potentials is not required to do this. Rather, only a few salient features need to 
be understood and adjusted appropriately. Beyond the obvious and well known ideas that the 
potential should be substantially, but not overly, asymmetric, and that the chromophore should 
absorb light well to the first excited state but only to some fraction of the sum-rule constraint, 
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these features are as yet poorly understood. 
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No. of Parameters Optimal f3 int h rnin /h m ax £io Aa^io 



1 


0.659525 









2 


0.701633 


1.19 x 10" 


- 1 




3 


0.701923 


2.25 x 10" 


-O 


0.554604 1.354420 


4 


0.708764 


3.45 x 10" 


r 
-0 


0.557740 1.268598 


5 


0.708928 


2.37 x 10" 


-7 


0.557752 1.269431 


6 


0.708836 


4.24 x 10" 


-6 


0.557762 1.269147 


7 


0.708951 


3.17 x 10" 


-7 


0.557763 1.269151 



Table I: For different numbers of parameters in the potential: (from the left) Optimized intrinsic hyper- 
polarizabilities for , together with the ratio of eigenvalues of the Hessian matrix of the intrinsic hyperpo- 
larizability. 



Eigenvalues of hessian 








Derivatives of dipol 


e matrix elements 


Numerically natural Measure 


Po 




PlO 






Ax' w 


-3.28 


-22.2 




-8.10 x 10" 


-i 


-3.21 x 10" 1 


2.15 


-3.14 x lO" 1 


-4.70 




-1.76 x 10" 


-i 


-7.07 x 10~ 2 


5.35 x 10~ 3 


-1.09 x 10~ 2 


-1.15 




-8.19 x 10" 


-2 


-2.05 x 10~ 3 


4.96 x 10~ 3 


-4.83 x 10" 3 


-2.47 x 10" 


-l 


-3.16 x 10" 


-2 


-1.47 x 10" 3 


1.13 x 10~ 2 


-2.58 x 10~ 4 


-8.03 x 10" 


-2 


-5.54 x 10" 


-3 


-9.80 x 10~ 5 


3.91 x 10" 4 


-1.66 x 10~ 5 


-4.65 x 10" 


-3 


-3.03 x 10" 


-3 


-1.23 x 10~ 5 


-7.71 x 10- 6 


-1.04 x 10~ 6 


-3.12 x 10" 


-4 


-2.02 x 10" 


-4 


3.44 x 10~ 7 


-4.85 x 10- 7 



Table II: Eigenvalues of Hessian for N = 7 parameters in each of three measures: (from the left) the 
numerically natural measure induced by the parametrization; the measure po = \tpo\ 2 and pio = \ipo\ 2 + 
\ipi\ 2 - Derivatives of the dipole matrix in the direction of the eigenvector associated with each eigenvalue 
are also displayed. 



16 




Figure 1: (Color online) (a) Ground and (b) first excited normalized wavefunctions corresponding to (c) 
optimized potentials for N = 2 to N = 7 parameters. Physically irrelevant positions and energies have 
been scaled out, consistent with eq. ( |18| ). The hyperpolarizabilities are noted in (a) and the ground and 
first excited energy levels are indicated in (c). 
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Figure 2: (Color online) Eigenvalues of the Hessian of /3i n t evaluated at its maximum plotted for different 
numbers of free parameters N. The dashed lines are intended as a visual aid only. 
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Figure 3: Variations of the potential and ground state wavefunction associated with eigenvectors of the 
Hessian matrix of f3i n t evaluated at the maximum; these have been adjusted so that each variation simulta- 
neously leaves E\—Eq = 1 and (x) = 0. The variations are plotted in descending order of eigenvalue. Each 
variation is calculated in two measures: the numerically natural measure (solid lines), and the measure of 
the ground state wavefunction (dashed lines). 



18 



